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Inspired by a recent discovery of a peculiar integer quantum Hall effect (QHE) in graphene, we 
study QHE on a honeycomb lattice in terms of the topological quantum number, with two-fold 
interests: First, how the zero-mass Dirac QHE around the center of the tight-binding band crosses 
over to the ordinary finite-mass fermion QHE around the band edges. Second, how the bulk QHE 
is related with the edge QHE for the entire spectrum including Dirac and ordinary behaviors. We 
find the following: (i) The zero-mass Dirac QHE (with a xy = :: f(2N + l)e 2 //i, N: integer) persists, 
surprisingly, up to the van Hove singularities, at which the ordinary fermion behavior abruptly 
takes over. Here a technique developed in the lattice gauge theory enabled us to calculate the 
behavior of the topological number over the entire spectrum. This result indicates a robustness 
of the topological quantum number, and should be observable if the chemical potential can be 
varied over a wide range in graphene. (ii) To see if the honeycomb lattice is singular in producing 
the anomalous QHE, we have systematically surveyed over square^->honeycomb <->7r-flux lattices, 
which is scanned by introducing a diagonal transfer t' . We find that the massless Dirac QHE 
[propto(2N + 1)] forms a critical line, that is, the presence of Dirac cones in the Brillouin zone is 
preserved by the inclusion of t' and the Dirac region sits side by side with ordinary one persists 
all through the transformation, (iii) We have compared the bulk QHE number obtained by an 
adiabatic continuity of the Chern number across the square<->honeycomb^ 7r-flux transformation 
and numerically obtained edge QHE number calculated from the whole energy spectra for sample 
with edges, which shows that the bulk QHE number coincides, as in ordinary lattices, with the edge 
QHE number throughout the lattice transformation. 
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I. INTRODUCTION 

Electrons on a honeycomb lattice, despite its simplic- 
ity, provide interesting problems in condensed matter 
physics, especially in its topological aspects. Field theo- 
retically, Dirac particles and associated gauge fields have 
been intensively investigated from a topological point of 
viewji^ so electronic properties for the honeycomb lattice 
may open new avenues for condensed matter phenomena. 
Indeed, there have been several proposals about peculiar 
properties in condensed matter systems that have zero- 
mass Dirac particles at low energy scales J»4 Apart from 
honeycomb lattice, zero-mass Dirac particles appear in 
condensed matter physics as effective theories in various 
guise. These include the d-wave superconductivity^, the 
so-called 7r-flux or chiral spin state^^ Anderson localiza- 
tion problems?^ spin related problems on the honey- 
comb lattice^S*ii and quantum phase transitions in two 
dimensions 

A seminal highlight, however, is the quantum Hall 
effect (QHE) in the honeycomb lattice which has 
recently been observed 16-17 in graphene, a monolayer 
graphite with a honeycomb array of carbon atoms. While 
the study of graphite has a long history, recent advances 
have been directed toward nanostructured graphite such 
as the carbon nanotube^ or nanographite with bound- 



ary magnetism expected to arise from edge statesii*!* 2 ^ 
In this context, the monolayer graphene is particularly 
interesting as an ideal realization of the honeycomb lat- 
tice, and the discovery of QHE has kicked off intenseive 
studies. 

Quite generally, topologically nontrivial states are 
characterized not by local order parameters as in sym- 
metry breaking states, but by geometrical phases; 21 ' 22 
where what are now known as topological orders can 
be realized i 2 ^* 2 ^^ 2 ^ One interesting consequence is that 
topological quantum numbers for the bulk can often be 
related with those for the edge states in finite systems. 
With this bulk-edge correspondence topological proper- 
ties which can be hidden in the bulk may thus become 
visible around the edges. A typical example is the edge 
states in QHE systemsi 27 ' 28 i 29 i 30 

For the QHE on the honeycomb lattice, we can then 
pose two fundamental questions: (i) While the low- 
energy theory around the band center (E = 0) is that 
of the zero-mass Dirac particle, which is now realized to 
give an anomalous QHE, how this should be taken over 
by ordinary theories as we go away from the band cen- 
ter? (ii) How should the bulk-edge correspondence look 
like for the zero-mass Dirac particle? The question (ii) is 
of special interest in the context of the zero mode of the 
massless Dirac particles^i*^ There, a bipartite structure 
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(chiral symmetry) in the honeycomb lattice is intimately 
related to the appearance of zero mode edge statesi^iiSi 

As for the question (i), Zheng and Andoi^ have nu- 
merically calculated the QHE on honeycomb with a self- 
consistent Born approximation. They have shown the 
anomalous QHE around E = 0, but they have also cal- 
culated the Hall conductivity for the entire energy region 
in this approximation. More recently, Sheng et al££ have 
examined the QHE in graphcne, computing the QHE 
number over the whole energy spectrum. They have 
shown that Dirac-like quantization of the Hall conductiv- 
ity appears only near the zero energy and conventional 
quantization can be observed in the band edge region. 
However, our question here is on the precise topological 
quantum number 38 in the region including around the 
boundary between the anomalous and ordinary ones. 

So in this paper, we explore the topological aspects for 
electrons on the honeycomb lattice in a magnetic field. 
We show the following: 

(i) The zero-mass Dirac particle behavior (with the 
Hall conductivity <j xy = ±(2N + l)e 2 /h, where N is an 
integer and we ignore spin degeneracy) persists, surpris- 
ingly, up to finite energies, at which the usual finite-mass 
fermion behavior abruptly takes over. The boundary en- 
ergies are identified to be the van Hove singularities. Here 
a technique developed in the lattice gauge theory enabled 
us to calculate the behavior of the topological number, 
which can become huge over the entire spectrum. This 
result indicates a robustness of the topological quantum 
number, and should be observable if the chemical poten- 
tial can be varied over a wide range in graphene. 

(ii) To see if the honeycomb lattice is singular in pro- 
ducing the anomalous QHE, we have systematically sur- 
veyed the systems that extend over square^honeycomb 
<->7r-flux lattices by introducing a diagonal transfer t' . 
We find that the Dirac region always exists (with its 
boundary dependent on t' and siting side by side with 
ordinary one) all through the transformation (except at 
the square lattice, which, with merged van Hove singu- 
larities, is rather singular in the present viewpoint). It 
implies the massless Dirac fermion forms a critical line in 
a parameter space. It does not occur as a critical point 
by adjusting parameters. Incidentally, we make an inter- 
esting observation, in the region honeycomb <->7r-fiux lat- 
tices, that the presence of multiple extrema in the band 
dispersion can give rise to a fermion doubling with the 
Chern number cx 2N rather than (27V + 1). 

(iii) We have then compared (a) the bulk QHE num- 
ber obtained by an adiabatic continuity of the Chern 
number across the square<-> honeycombs 7r-flux transfor- 
mation and (b) numerically obtained edge QHE number 
calculated from the whole energy spectra for sample with 
edges. The result shows that the bulk QHE number coin- 
cides, as in ordinary lattices, with the edge QHE number 
throughout the lattice transformation. Incidentally, the 
E = flat edge mode persists in strong magnetic fields. 

The organization of the paper is as follows: In the 
next section, we define the tight-binding model on the 



honeycomb lattice, where we introduce a diagonal trans- 
fer to tune the position of the van Hove singularities to 
go over to the square and 7r-flux lattices. We present 
numerical results for the energy spectra as a function of 
a magnetic field (Hofstadter diagram), which enables us 
to infer a topological relationship among the honeycomb, 
square and 7r-flux lattices. In Sec. IIIII we compute the 
Hall conductivity of the bulk as a function of the chemi- 
cal potential, based on a lattice-gauge theoretical method 
developed in Ref4?-. This calculation is manifestly gauge 
invariant and guarantees integer Chern numbers, so that 
the method is powerful in evaluating the QHE topological 
number over the whole spectrum, including the van Hove 
singularities, which turn out to accompany singular be- 
haviors in QHE. There the topological equivalence in Sec. 
Ill Bl is thus confirmed with respect to the bulk topologi- 
cal properties. In Sec. IIV1 we show that the conversion 
square<->honeycomb <->7r-flux has a virtue of enabling us 
to derive the Diophantine equation for the Dirac-fermion 
regime from the adiabatic principle for the topological 
quantum number, which exists as far as the energy gap 
in which Ep reside does not close. Sec. [V]is devoted to 
the edge states of the model. Edge states in finite (cylin- 
drical) systems are analyzed with the transfer matrix for- 
malism, and we show that edge states, with Laughlin's 
argument^ indeed reveal the Dirac like quantization sit- 
ting next to the conventional quantization separated by 
van Hove singularities. Thus, the bulk-edge correspon- 
dence is confirmed. A summary and discussion is given 
in Sec. EH 



MODEL AND TOPOLOGICAL 
EQUIVALENCE 



II. 



Let us start with defining the Hamiltonian for a model 
that interpolates honeycomb lattice with square and tt- 
flux lattices. This parameter plays a key role in exam- 
ining the topological properties in terms of the adiabatic 
principle. 



A. Model 

The Hamiltonian for two-dimensional tight-binding 
systems in a uniform magnetic field B applied normal 
to the plane is 

E -^honeycomb ~t~ E . 

Here, -Hhoneycomb is the tight-binding model on the hon- 
eycomb lattice with nearest-neighbor hopping, 

^honeycomb =t £ [4 (j)c (J) + 6°"+* c\ (j)c (j - e 2 ) 



-cl(j + ei)c (j) 



h.c, 



where i = — 1 is the transfer (taken to be the unit of 
energy). Since the honeycomb, a non-Bravais lattice with 
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two sites per unit cell is bipartite, we have defined two 
kinds of fermion operators, c (j) and c. (J) as in Fig. 
^ where j = j\e\ + fi e 2 labels the unit cell with two 
translation vectors e\ = (3/2, \/3/2)a and e 2 = (0, \fZ)a. 
The magnetic field is characterized by 

cj> = BS 6 /(2ir)=p/q, 

the magnetic flux (assumed to be rational) penetrating 
each hexagon of area S 6 = {3V3/2)a 2 . we have adopted 
the Landau gauge for the vector potential. For simplicity 
we neglect the spin degrees of freedom, since the graphene 
has a very small Zeeman splitting. If we add a hopping 




FIG. 1: A honeycomb lattice with a unit cell and a extra 
transfer t' indicated. 

across each hexagon, 

H' = t'Y,e- l2 ^ Ul+1/2) cl(j + ei- e 2 )c (j) +h.c, 



we shall stress, it is between the two singularities that 
a zero-mass Dirac like gapless dispersion appears in the 
honeycomb lattice. We can in fact rigorously showSi that 
the zero gap with a linearly vanishing density of states 
around E = is not an accident for honeycomb, but 
exists for — 3 < t' < 1, so such a situation persists down 
to the 7r-flux lattice. 

Ordinary band dispersions and density of states reside 
outside the van Hove singularities. So our question here 
is how the quantum Hall effect should look like along the 
sequence, square<-^honeycomb <-+7r-flux. 



B. Topological Equivalences 

In a magnetic field B, the spectra of the lattice Hamil- 
tonians against B oc <f> take fractal shapes, usually called 
Hofstadter's diagram^ where hierarchical series of en- 
ergy gaps exist. The butterfly is deformed as t' is var- 
ied. Here, an adiabatic principle plays a crucial role: it 
dictates that one can keep track of quantum mechani- 
cal ground states when the Hamiltonian is transformed 
with a continuous change of parameter(s) if there is a 
gap above the ground state and if the gap remains finite 
throughout*^ Let us apply this argument to the present 
model, where t' is the adiabatic parameter. In Hofs- 
tadter's butterfly we can see many Landau bands. The 
change of t' makes some of the bands merge (or, more 
precisely, some of the gaps between Landau bands merge, 
since the spectrum is fractal). The Landau levels are 
characterized by the quantum Hall numbers, which are 
topological (Chern) numbers as will be discussed in Sec. 
IIII and the numbers remain unchanged against the adi- 
abatic change as long as the gap in which Ep reside does 
not collapse. This is a topological stability. 



we can change the system continuously from t'/t = — 1 
(which is referred to as 7r-flux lattice) to (honeycomb) 
and 1 (square). 

In the momentum space, the Hamiltonian for cf> = p/q 
is expressed as 

where c^(fc) is a 2q dimenisional vector and h(k) is a 
2q x 2q matrix (see Appendix |A"|) . In zero magnetic field, 
we have 

h(k)-tf ° A(fe) 

A(fc) =1 + e lk2 + e lkl [1 + (t'/t)e- lk2 \ . 

To see how t 1 shifts the position of van Hove singularities, 
we show the energy dispersions and the density of states 
in the zero magnetic field in Figs. [2] and [31 The square 
lattice has a van Hove singularity at the band center, 
which splits into two as we decrease t'/t from 1. As 



1. Topological equivalences between van-Hove singularities 

The Hofstadter diagram for honeycomb lattice was first 
obtained in Reff2S. Rammal has already noted the pres- 
ence of the E = Landau level which is outside Onsager's 
semiclassical quantization scheme. The spectrum for the 
honeycomb lattice is symmetric about (f> = tt. In Figs. 

we show how the spectrum versus <p oc B is deformed 
with t' for the square lattice (t'/t — > 1), or for the ir- 
flux lattice (t'/t — > —1). In the present construction the 
flux per plaquette on the square lattice corresponds to 
half the hexagon in the honeycomb lattice, which is why 
Hofstadter's diagram has a period of 4ir for t' ^ 0. The 
property that the spectrum is invariant against (j) — > — 4> 
as well as against t' <-> — t', 4> «-> 1 — is retained through- 
out. The honeycomb lattice with t' = corresponds to 
the self-dual point for the t' <-> — t' symmetry. 

We can immediately notice from Fig. Eta)(b) that the 
large gaps just above and below the zero energy remain. 
This implies that the topological number should remain 
the same when E-p lies in the gap. A closer examination 
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(a) x flyx: (■/(=-■} [h] t ft=-Q.5 (c) hnn#y$nmb: f/t=Q (d) f$=+Q.$ (r) square : f:t=+ 1 




FIG. 2: Energy dispersions for (a) t'/t = -1 (yr-fiux lattice), (b) t'/t = 0.5, (c) t'/t = (honeycomb), (d) t'/t = 0.5, and (e) 
t'/t = 1 (square). 



(a) Kflux:t'/t=-1 (b) f'/f=-0.5 (c) honeycomb: t'/t=0 (d) f'/t=+0,5 (e) square : t'/t=+1 




■4 -2 2 4 -4 -2 2 4 ^4 i 2 4 -4 -2 2 4 -4 -2 2 4 



FIG. 3: Density of states for (a) t'/t = -1 (tt-Aux lattice), (b) t'/t = 0.5, (c) t'/t = (honeycomb), (d) t'/t = 0.5, and (e) 
t'/t — 1 (square). 



shows, surprisingly, that this holds for other gaps, all 
the way up to a finite energy, E c , in fact. Figure |S| is a 
blowup of the low magnetic field region. We can see that 
there is no level crossing for other energy gaps as well 
in the energy region E < E c , so that we can keep track 
of the main gaps in the adiabatic change between the 
honeycomb and 7r-flux (t'/t = <-> —1). The boundary 
energy E c is seen to reside around the region where the 
Landau level fan extending from E = meets another 
Landau fan extending from the band edges \E\ ~ 3. We 
shall show below, by using the calculation of the quantum 
Hall numbers, that we can identify E c as the van Hove 
energies (which are functions of t'). 



2. Topological equivalences near band edges 

The above situation is exactly symmetric for the 
crossover of the honeycomb <-> square (t'/t — <-> 1), 
where we have only to replace (f> with 1 — 0. This means 
that, in the weak magnetic field region depicted in Fig. 
large gaps around E = are closed as we go to the square 
lattice, while we can keep track of the Landau fan start- 
ing from the band edge. This property is in fact highly 
nontrivial, since the Landau levels is usually analyzed in 
the effective-mass approach only near the band center or 
band edges. This is exactly why we want to confirm the 
topological equivalences in terms of the quantum Hall 
number in the next section. 



III. TOPOLOGICAL EQUIVALENCES IN THE 
BULK 

A. Hall condactance 

It is now firmly established that the Hall conductivity 
of the noninteracting 2D electron systems, as described 
with the linear response theory^ may be regarded as 
a topological quantum number42*2i Namely, when the 
Fermi energy lies in the n-th gap, the Hall conductivity 
is given by 

e 2 

o~ xy = -~^c F (E F ), 

where cf is an integer called Chern number, which de- 
scribes how the wave function (a vector bundle) responds 
to a vector potential generated by Berry's gauge poten- 
tial in the Brillouin zone. To compute the number, we 
need to diagonalize the Hamiltonian 

h(k)<pj{k) = e,-(fe)^(fe), 0' = 1.2,- • • ,2q), 

where h(k) is the 2q x 2q matrix defined in Eq. ljAip . 
<Pj(k) an eigenvector, and we assume that t\(k) < • • • < 
e„(fc) < Ep < e n+ i(k) < ••• < £2q(k) holds over the 
entire Brillouin zone. We can then define Berry's gauge 
potential, 

A ij (k) = <pl(k)V k < Pj (k), (l<i,j<n), (2) 

where Ay(fe) is an anti-Hermitian n x n matrix. The 
Chern number is given as the U(l) part of the above 



FIG. 4: Hofstadter's diagram (energy spectrum against the magnetic flux (j>) for (a) t'/t — —1 (7r-flux lattice), (b) t'/t — 0.5, 
(c) t'/t = (honeycomb), (d) t'/t = 0.5, and (e) t'/t = 1 (square). 




FIG. 5: Blowup of Hofstadter's diagram in a weak magnetic field region for (a) t' /t = —1 (7r-flux), (b) t'/t = 0.5, (c) t'/t — 
(honeycomb), (d) t'/t = 0.5 and t'/t = 1 (square). The arrows indicate the positions of van- Hove singularities. The blue lines 
indicate positions of a flux (f> — 1/31, which corresponds to the one in the Figs© 



U(n) gauge potential, 



c F (E F ) = — / Tr dA, (3) 
2m J 



where A(k) = A fl (k)dk f2 is a one-form . 

This formulation, a non-Abelian extension of Berry's 
gauge potential's*^ may seem too elaborate, but is use- 
ful when there are multiple bands below the Fermi en- 
ergy. Namely, this formula holds even if some gaps in the 
Fermi sea are closed, as long as the bands in question, 
e„(fc) and e I1+ i(fc), do not cross. In the special case when 
all the bands in the Fermi sea are separated with each 
other, the Chern number, Eq.Q, is simply the sum of 
the Chern numbers assigned to individual bands^i but 
we have opted for the above formula, since some of the 
gaps do collapse as t' is varied, as we shall see. 



B. Lattice gauge theory technique 

To actually evaluate the topological integer for the 
present system, we need to calculate eigenfunctions of 
the 2q x 2q Hamiltonian (|A1|) . Diagonalization can be 
done only numerically, so that the eigenfunctions are ob- 
tained in practice on mesh points in the Brillouin zone. 
This can cause a serious problem if we want to obtain 
the Chern numbers, especially around the van Hove sin- 
gularities, since they can behave wildly there. The prob- 
lem becomes even more formidable if we want to main- 
tain the calculation manifestly gauge invariant and pin 
point the Chern numbers as integers. Here we adopt a 
method^iMi developed recently in the context of the lat- 
tice gauge theorym^^^^^^ This implementation pre- 
cisely guarantees manifest gauge invariance and integer 
Chern numbers. 

Let us first compute the eigenfunction of the Hamilto- 
nian on meshes in the Brillouin zone. In what follows, 
we choose the Landau gauge for the magnetic field, so 
that the mesh points are denoted as (kj t , kj 2 ), where 
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FIG. 6: The Chern number (Hall conductivity in unit of —e 2 /h) for magnetic field <f> — 1/31 is plotted against the Fermi 
energy E F : for (a) t'/t = -1 (yr-flux), (b) t'/t = 0.5, (c) t'/t = (honeycomb), (d) t'/t = 0.5, and (e) t'/t = 1 (square). The 
red lines indicate the steps of two in the Chern number sequence (a xy — ±(2N + l)e 2 /h, N: integer), while the blue lines step 
of one (a xy = ±Ne 2 /h). 



discrete sets of momenta with — 0, . . . , — 1 are 
denned by kj 1 — 2irji/(qNi) and kj 2 — 2^32/ N 2 in 
the Brillouin zone extending over < k\ < 2n/q and 
< lc2 < 2ir. Diagonalization of h(k) on these sites 
k = kg = (kj 1 ,kj 2 ),l — 1, . . . , NiN 2 provides eigen- 
functions with h(kg)ipj(k() — ej(ki)(pj(kt), which can 
be chosen to satisfy the periodic boundary condition 
tpj(ke + N^ft) = ifj{ki), where p, is a vector of length 
|i| = 2Ti/{qN 1 ) and |2| = 2tt/N 2 along fc M . 

Provided that the Fermi energy lies between the nth 
gap (namely, e n (ke) < E F < e n+ i(k e ) for every k e ), and 
hence the Fermi sea is composed of n bands (some of 
which may merge), we define a U(l) link variable of the 
Fermi sea as 

U^kt) = IdetU^kt)]- 1 debU^kt), 

where 

(U^ij = tptik^tpjiki + A), (1 <i,j< n). 

The link variables are well-defined except at singular 
points det U^kg) = 0, which corresponding to "vortices" 
and "antivortices" in the wave function, and we can al- 
ways make the mesh avoid them with an infinitesimal 
shift. With the link variable we can define a lattice field 
strength by 

F 12 (k e ) =lnC/ 1 (^)[/ 2 (fef + i)C/i(fe^ + 2)- 1 [/ 2 (^)- 1 , 

where the principal branch of the logarithm with — 7r < 
Fi2(kg)/i < 7r is taken. The field strength is by definition 
invariant under gauge transformations, so that no specific 
gauge fixing is required. The Chern number on the lattice 
is defined as 

c F (E F )~^F 12 (kz). (4) 

t 

The Chern number thus computed is strictly an integer, 
since F±2 just counts the number of vortices minus 



the number of antivortices in the Brillouin zone (i.e., the 
number of times the branches in In are traversed) , as has 
been proved in Ref4£. Therefore, it should give the num- 
ber in the continuum, cp appearing in Eq. @ when the 
mesh in the Brillouin zone is sufficiently dense. This is 
due to the integral character of the lattice Chern number: 
For large N^, we have 

U^k) = l + Tr^(fe^fc M , 

with defined by Eq.©, and we have 

F12 ^InHil + TrA^ + OdSkl 3 ) 
□ 

=Tr [d x A2{k) - d 2 A 1 (k)} SkiSk 2 + 0(\Sk\ 3 ). 

When coupled with Eq. Q, this reduces to Eq. |(3}, the 
Chern number in the continuum Brillouin zone. 



C. Dirac vs ordinary fermion quantization 

Let us now show the numerical results for the Chern 
number cp in Eq. Q as a function of the chemical po- 
tential E F for various values of t' . While the integral 
Chern number c F is defined only for each gap, we have 
plotted cf as step functions, which makes sense as long 
as the magnetic field is not too large, as in Fig. Elwith 
4> = 1/31, since the width of each Landau band is then 
much smaller than the size of gaps. We have displayed 
the values for spinlcss fcrmions: 

The result shows a striking answer to one of our key 
questions: what is the fate of the Dirac-like behavior as 
we go away from E = 0. When the Fermi energy is 
swept in the honeycomb lattice, the Dirac-like Hall con- 
ductivity steps of two (or four when spin degeneracy is 
included) in units of —e 2 /h exist around E = as has 
been noted by many authors>^*iSiiL^i Let us call this 
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Dirac-Landau behavior. Now, we can see that this Dirac- 
Landau behavior persists all the way up to some energy, 
which we identify to be the van Hove singularities ap- 
pearing in Fig. |3J At these energies we have then a huge 
step accompanied by a sign change in the Hall conduc- 
tivity. This result implies the following: The effective 
theory near the zero energy is Dirac-like fermions in the 
continuum limit, for which an unconventional quantiza- 
tion of the Hall conductivity has been derived from the 
Dirac Landau levels The present calculation leads 
to the conclusion that the unusual property extends to 
an unexpectedly wide region of energy in the actual lat- 
tice fermion system. Outside the van Hove energies (i.e., 
in the band-edge regions), we recover the conventional 
QHE where the step changes by one in units of — e 2 //i 
(which we now call Fermi-Landau). This implies that a 
huge step accompanied by a sign change has to occur at 
the boundary between Dirac-Landau and Fermi-Landau 



behaviors, or the bands just at the van-Hove singularities 
induce a change, which is topological in that the relevant 
quantum numbers are topological. 

A second striking feature in Fig. [fj]is that the Dirac- 
like Hall conductivity oc (2 A + 1) persists all along 
— 1 < t'/t < 1 in the 7r-fluxs honeycomb<->square se- 
quence, except at the square lattice. Namely, there al- 
ways exists the Dirac-Landau behavior between the (in- 
nermost) van Hove singularities, except for the case when 
the singularities merge at t' — 0. We can check from 
the band dispersion (Fig|2J that in this whole region we 
have two zero-mass Dirac cones in the Brillouin zone, 
which should cause this persistence of the anomalous 
QHE number. 

We can now summarize the Hall conductivity for the 
spinless fermions on the honeycomb lattice for the entire 
energy region as 



he 



-(N+l) E F <-|t|, AT = 0,1, 2, 

-(2JV+1) -\t\<E F <0, A = 0,1, 2, 

-(2JV+1) 0<E F <+\t\, A = 0,1, 2, 

■{N + l) \t\ < E F , A = 0,1, 2, 



Landau level index counted from the band bottom 
Dirac-Landau level index counted from the zero energy 
Dirac-Landau level index counted from the zero energy 
Landau level index counted from the band top 



We should again double these numbers if we include the 
spin degeneracy. While the unconventional quantization 
around the zero energy has been beautifully observed 
experimentally 1617 , pushing the chemical potential fur- 
ther to approach the van Hove energies should detect the 
topological phase transition. 



We can note that, to be precise, there are in fact 
two sets of van Hove singularities (i.e., four in total, see 
FigEJh)) in the honeycombs 7r-flux region. There, we 
have the Dirac-Landau behavior in the region between 
the innermost singularities, Fermi-Landau behavior in 
the next region, and another behavior outside the outer- 
most singularities with the doubled Chern number oc 2A. 
Although this may seem a trivial detail, we can in fact 
make an interesting observation. Namely, we can raise 
a question: how the band-edge region can have doubled 
Chern number than Fermi-Landau behavior? We can 
identify the origin of the doubled Chern numbers in the 
fermion band as the presence of two minima (maxima) in 
the dispersion in the Brillouin zone at the band bottom 
(top). So we can make an observation that a degenerate 
minima (or maxima) in the band dispersion can give rise 
to a Fermion species multiplication. We shall reinforce 
this view in terms of the topological equivalences in Sec. 
Ill Bl and edge state pictures in Sec. [V] 



D. Degeneracy of Landau levels and the lattice 
effect 

Let us now look more closely at the distinction between 
the Dirac-Landau and Fermi-Landau levels in terms of 
the degeneracy of the Landau bands. The physical reason 
why we have a step of two in the quantized Hall conduc- 
tivity for the Dirac-Landau levels can be ascribed, as in 
many other literatures, to the degeneracy (fermion dou- 
bling) of the bands. Let us take an example of <p = 1/31 
employed in Fig. EJc), for which we have 2q = 62 eigen- 
values tj{k). Among them, there are 7 pairs of twofold 
degenerate levels around E = with a total Chern num- 
ber of 2, while other 48 levels are nondegenerate. Two of 
the nondegenerate Landau levels are located at the van 
Hove energies and carry a huge Chern number of —30. It 
is such large numbers that require the accurate estimate 
with the present, lattice-gauge numerical algorithm. 

The level adjacent to these have the ordinary Chern 
number of 1 each. So this abrupt change in the topolog- 
ical number defines the clear boundary. To be precise, 
however, we can note that approximate degeneracy in 
each two-fold Dirac-Landau levels (which is exact for a 
Dirac fermion with a rigorously linear dispersion) become 
slightly lifted toward the van Hove energy, while the gap 
across the next two-fold level becomes smaller, where the 
splitting (0.05t for <fr = l/31)and the gap (0.02i) become 
comparable (which is still much smaller than the thick- 
ness of the line that plots the energy dependence of the 
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Chern number in FigHO so the figure still makes sense). 
The deviation becomes more pronounced for stronger 
magnetic fields. We shall discuss this feature in Sec. E] 
(see, e.g., Fig. [TJ3 (a) for <j) — 1/5), where the gap can 
become comparable with the Landau band width. 



IV. DIOPHANTINE EQUATION 

Having numerically calculated the Chern numbers, we 
can now raise the question: can we calculate these num- 
bers algebraically? For square lattice systems, the Dio- 
phantine equation according to Thouless et al^ is known 
to be a simple yet powerful tool to compute the quan- 
tum Hall number of the Fermi seas. For the honeycomb 
lattice, however, the Diophantinc equation has not been 
obtained. This is exactly where we can exploit the topo- 
logical equivalence establised in the previous sections to 
calculate the Chern numbers by the use of the Diophan- 
tine equation. Natural interest here is how the Dirac- 
Landau quantum Hall number of (27V + 1) comes about. 
So we start with the Diophantine equation for the square 
lattice, where the Chern number cj of the Fermi sea com- 
posed of J bands with flux per plaquette $ = P/Q is 
given by the formula 

J=Pcj (mod Q), \cj\ < Q/2. (5) 

Now let us make its adiabatic continuation to the hon- 
eycomb lattice with magnetic field <p = !/?■ Since the 
7r-flux lattice is constructed here from the honeycomb 
lattice by the introduction of t'/t = — 1, the flux per 
square plaquette corresponds to (f>/2 for the hexagonal 
plaquette with the shift in the period discussed in Sec. 
Ill Bl Therefore, total flux per square plaquette of the 
7r-flux lattice, $, reads 

_ P 1 <j> q + l 

$ = — = h — = 

Q 2 2 2q ' 

where the 1/2 flux in the above equation is due to the 
7r-flux. Let us first take the case of q+l: prime (q: even). 
Then we have 

P=q + 1, Q = 2q. (6) 

Equation (J5J together with Eq. © is a Diophantine 
equation for the honeycomb lattice with even q valid in 
the Dirac fermion regime. To identify J we assume for the 
honeycomb lattice that there exist a zero energy Dirac- 
Landau level, and the Dirac-Landau levels between the 
van Hove singularities are twofold degenerate. 

Now let us calculate the Hall conductivity around the 
zero energy for the honeycomb lattice when E-p lies in the 
(iV + l)-th gap from the center (e.g., N — is the first 
gap above the zero energy). Note that the total number 
of the bands is 2q, so that the number of the negative- 
energy bands is q — 1. Hence the number of bands below 



the (iV+l)-th gap is q-l + 2(N+l), i.e, J = q + 2N + l. 
From Eq. (JSJ, the Diophantine equation then reads 

q + 2N +l=(q + l)cp, (mod 2q), 

or 

(g+l)(cp-l)=2JV = 2JV(g + l), (mod 2q). 

This has an obvious solution, 

c F =2N + 1 : E F > 0, 

and its particle-hole transformed cf = —{2N +1) for 
E-p < 0. Remarkably, this precisely gives the Chern num- 
ber of the Dirac fermion, as computed in the previous 
section. 

For an odd q, we put 

q =2n + 1, n = 1,2, • • • . 

Then the mapping to the 7r-flux lattice reads 

_ P _ q+l _ n + 1 
~ Q 2q ~ 2n + 1 ' 

Namely, we have 

P = n + 1, Q = 2n+1. (7) 

Note that the total number of the bands in this case 
is Q = q = 2n + 1 rather than 2q in the adiabatically 
transformed tt flux system. Then, the number of bands 
in the negative energy is n, and the number of bands 
below the (N + l)-th gap is n + (N + 1), i.e, 

J =n+l + N, 

which gives, via Eqs.JSJ),©, a Diophantine equation, 

n + l + N =(n + l)cp, (mod2n + l), 

or (n + l)(cp — 1) = N. This expression is equivalent to 
— n(cp — 1) = —2nN. This has a solution cp = 2N + 1, 
which gives the same Chern number for the prime q + l 
case. 



V. EDGE VS BULK FOR THE HALL 
CONDUCTIVITY IN THE HONEYCOMB 
LATTICE 

When a system has a nontrivial topological structure, 
they should generically exhibit characteristic properties 
related to edge states. For the honeycomb lattice, edge 
states and their flat dispersion have been intensively dis- 
cussed with its relevance to spin alignmen1ji2iS£. Another 
way of saying is that the honeycomb system has a bi- 
partite symmetry, which guarantees the presence of the 
zero mode edge states^, i.e., there exist macroscopic edge 
states at the zero energy unless the bipartite symmetry 
is broken. This gives rise to dispersionless edge states. 
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The flat bands are unstable against bipartite symme- 
try breaking (Peierls instability), which is, e.g., realized 
when an antiferromagnetic spin ordering occurs. This 
may be viewed as a topological origin of the boundary 
spin moments. 36 

Now, the problem of edge states vs bulk states is par- 
ticularly interesting for the QHE, since the problem of 
how the Dirac-Landau QHE oc (27V + 1) is related to 
the edge states is fundamentally interesting. One of the 
present authors has shown, for 2D square (or anisotropic 
square) systems, that the edge states whose energy dis- 
persions run across the Landau bands have topological 
QHE numbers, for which 

edge _ -bulk 
xy xy 

by identifying the connection between the topological in- 
tegers for the bulk and for the edge statesi 3 ^^ So let us 
now focus on the edge states of the honeycomb lattice in 
terms of the topological structure. It also allows us to 
clarify the bulk-edge correspondence. 



A. Transfer matrix formalism 



(a) 





1. Transfer matrix 

We follow the analysis in re fs22iiffii£i to reduce the sys- 
tem to a one-dimensional model by making a partial 
Fourier transform of the fcrmion operators in 2-direction, 



Co,, (j) 



1 



c ,.{ji,k 2 ). 



(8) 



This yields a ^-dependent series of one-dimensional 
Hamiltonian 

Hm(k 2 ) = ^o.Cn,fe)4(ji, fc 2)c.(ji,fc 2 ) 



+ £.o(jl,fc2)4(j'l + l,fc 2 )c (jl, fa) 



h.C. 



(9) 



where fe-dependent hopping parameters are 
to.(h,k 2 )=t(l + e ik2 - i2 ^), 

t.a{jl,k 2 ) =t [l + {t'/t) e ifc2-i27T0(i 1+ l/2) 



Expanding one-particle eigenstates \E, k 2 ) with energy E 
and momentum k 2 as 



\E,fa)=J2 



^.(S,j 1 ,fe 2 ) C t(i 1 ,fe 2 )|0) 
+ MEJuk 2 )ct(j ll k 2 )\0) 



we find that the Schrodinger equation Hi^,{k 2 )\E, k 2 ) 
E\E,k 2 ) is cast into a matrix form, 



V'o(jl) 

V'.O'i) 

^o(ii) 



Mo.(?i 
-- M.o(ji) 



^o(7l-l) 



with 



Therefore, we have 



E 


t.oO'l-l) 




tS.O'i) 


l 





E 


to.(jl) \ 


KoUi) 




1 


) 



(10) 



tpo(jl) 



M t {j x ) 



FIG. 7: A cylindrical system with zigzag and bearded edges. with 



Mt(ji) = M.o(ji)M .(ji 



(11) 



where every quantity is a function of (E,ji,k 2 ). They 
are a set of equations for 2L\ — 1 variables t^»(j x )(jx = 
1, • • • , L x - l),ipo(jx)(jx = 1, • ' ' ,L X ) with open bound- 
ary conditions tp o (0) = ip {Li) — (See FigQl. To im- 
pose this condition, we consider cylindrical systems with 
a zigzag edge at one end and a bearded (or Klein's) edge 
at the other as illustrated in Fig. This is just a techni- 
cal convention to apply the following transfer-matrix for- 
malism. In terms of the energy spectrum, zigzag-bearded 
system has an E — edge state in zero magnetic field 
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FIG. 8: The energy spectra for the cylindrical system in Fig. Qwith the flux (f> — 1/5 for (a) t'/t = —1 (7r-flux), (b) t'/t = 0.625, 
(c) t'/t — (honeycomb), (d) t'/t = 0.5, and (e) t'/t — 1 (square). Blue (red) lines are edge states localized at the bearded 
(zigzag) edges, while the bulk energy bands are shown as grey regions. The topological number, I(E-p), when E-p is in a major 
gap in the negative energies are: (a) 2, — 1; (b) 1, 2, 3, —1; (c) 1, 2, 3, —1 (d) 1, 2, 3, (ill-defined) and (e) 1, 2, 3, 4. A 
topological transition accompanied with a discrete Chern number change by Ac = —5 (from —4 to +1) at the 4-th major gap 
occurs between t'/t — 0.5 — > 1. The green lines in the right are the topological numbers, I(Ep), obtained by counting the edge 
states. 



that extend over the entire Brillouin zone, as was first 
pointed out In passing, we mention that the local 

transformation 

V>. (E, ii) -> - i). (E, h ) = ip.(-E,h), 

ME,Ji) -^t>o{E,h) = M~E, 3l ), (12) 

yields the same Schrodinger equation but with — E. This 
is due to the bipartite symmetry of honeycomb: Eigen- 
states for ±E are paired with the same amplitudes 
\i/j»,o(jn)\ except for the zero energy. 

Here note the periodicity M s (ji + q) = M s (ji) for s = 
o», s = »o, and hence M t (j\ + q) = M t (ji). We can now 
define the transfer matrix for a unit period, 

M(E,k 2 )= Y[M t {E, Ju k 2 ), (13) 

then 

(*■&>" )-"($$)■ ( "> 

For L\ — ql with an integer i we have 

\ MLi) ) W°(°) )' 

with the boundary condition ip o (0) = tpo(Li) = and 
•0.(1) = 1, which is equivalent to an algebraic equation 
(M )i\{E) = for E, which is a polynomial of order 
2L± — 1 and has as many real roots. 

We can extract the edge states from the whole 
spectrum 29 ! 30 . Namely, 2q — 1 roots of M2x{E) = 0, a 
polynomial equation of order 2q — 1 , give the edge state 



energies, Ej{ki) for j = 1, • • • , 2q — 1. By contrast, usual 
cylindrical systems where both edges are zigzag, bearded, 
or armchair do not allow such extraction. For exam- 
ple, in the case for zigzag edges, the boundary condition 
is replaced by tp o {0) — tp,(Li + 1) = in the present 
formulation. Hence, the total spectra is determined by 
(M e )n(E) = 0, giving 2L\ real roots. Edge states are 
contained in these spectra, but there is no simple way to 
directly extract them only. 

A remarkable feature here is that, on top of the edge 
states across adjacent Landau bands, there always exist 
exactly constant E — edge states, since the number 
of the edge states are odd 2q — 1. Namely, a system 
with zigzag edges has a zero-energy edge mode in some 
region of k%, a system with bearded edges in the other 
region, and the present system has the mode over the 
entire region. This is a simple extension of the discussion 
in the absence of a magnetic fieldi^ 3 - 

We can further identify the spatial position of edge 
states: With Mu(Ej), we can show that 

\Mu(Ej)\ < 1 : localized at the zigzag edge, 
|Mii (-Ej-)| > 1 : localized at the bearded edge. 

With Laughlin's argument^, supplemented by the above 
behavior of the edge state (as a function of the momen- 
tum £2)5 one can identify a topological number I (E-p), 
the number of electrons carried by Laughlin's adiabatic 
procedure, for each edge state branch when the chemical 
potential E-p traverses the branchiSLSS^SiSi 

Having determined the spectra for the edge states with 
the transfer matrix, we can now relate the edge states to 
the Bloch functions of the bulk systems: Let us relax the 
open boundary condition, and impose instead a periodic 
boundary condition H(j±-\-q, k%) = H{j\, k 2 ) on the local 
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Hamiltonian Then the Bloch theorem for the bulk 
state leads to + l) = e tkiq ip»(ji) an d V>o(ji + q) = 
e ikiq ip (ji), or 

/ ^.(il + q + 1) \ Jfcl9 / ^.(ii + 1) \ 

*v ^o(ii+g) J ^ Mn) )' 

Comparing this with Eq. 114(1 , we see that the eigenvalue 
of the transfer matrix M is e lkiq , i.e, just a phase factor. 
Therefore, the energy bands are specified by the following 
condition for the eigenvalue p of the transfer matrix M , 
det(/ol - M) = 0, with \p\ = 1. As shown in the App. [Bj 
the condition for the energy band is given by |TrM| < 2. 




FIG. 9: A Riemann surface E2 q _i(£;2) f° r the Bloch function 
which is a complex energy surface for defining p(E) consis- 
tently. Its genus is 2q — 1, which coincides with the number 
of energy gaps. 

The eigenvalue of the transfer matrix is also written as 
p—\ \TtM - [(TrM) 2 - 4detM] 1/2 |, where we have 
to fix the branch of the square root. Then the Riemann 
surface, £29- 1, of the Bloch function is given by that of 
the eigenvalue p(E) and its genus is generically 2q — 1 
which are the number of energy gaps (FigEJ. With this 
Riemann surface and the arguments in RefsSLiffiiSi, one 
can establish the bulk-edge correspondence as 

c F (E F ) =I(E F ), 

where c F (Ep) is the Chern number obtained in section 
III, while I(E F ) is the topological number for edge states 
which corresponds to the number of intersections on the 
Riemann surface X 29 _i. 

B. Edge-bulk equivalence and topological 
transition 

Before we discuss the experimental situation, let us 
look at the case of strong magnetic fields with <f> = 1/q 
with q = O(l). This case, while unrealistic, is useful for 
heuristic purposes. In Fig. 00 the spectra of the cylindri- 
cal systems are shown for a relatively large <f> = 1/5. The 
shaded regions are Landau bands, while the blue (red) 
curves across the bands correspond to edge states local- 
ized at the bearded (zigzag) edges. It should be noted 
that we always find a dispersionless edge state at the zero 
energy, which is due to bipartite symmetry. This is the 



state employed for the boundary magnetic moment, but 
does not contribute to the topological number, since the 
number is ill-defined due to the gap-closing at the zero 
energy. 

For the honeycomb lattice there are 9 bands for <fi — 
1/5, as observed in Fig. 03(c). By counting the number of 
the edge states, we can see that the topological numbers 
I(E F ) are 1, 2, 3, -1, +1, -3, -2 and -1, when E F lies 
in the jth major gap for j = 1, • • ■ ,8, respectively. Note 
that gap-closing occurs at the zero energy, and the incre- 
ment of the topological number across these degenerate 
bands is AI = 2 which corresponds to the Dirac-Landau 
level quantization of the Hall conductivity discussed in 
Sec. HTU 

If we move on to 7r-flux lattice (t'/t : — > —1), the 
Dirac fermion level with AI = 2 at the zero energy per- 
sists (Fig. 00(a) and (b)). This is the edge-bulk correspon- 
dence version of the topological equivalence (I) in Sec. 
Ill Bl If we move on to square lattice (t'/t : — > 1), on the 
other hand, we can observe a topological change in these 
bands: At t'/t — 1, the topological numbers I(E F ) reads 
1, 2, 3, 4, —4, —3, —2, — 1 versus E F , which dictates that 
E = is now a van-Hove energy with | A/| = 8 ^> 1. This 
corresponds to the fact that van-Hove energies merge to 
wash out the Dirac-Landau region precisely at t'/t = 1, 
so the square lattice is singular in this respect. The merg- 
ing of the relevant Landau bands (4th-6th for = 1/5) is 
clearly seen in the spectra at t'/t = 0.625 shown in Fig. 
|Sfd), where the topological number —5 moves to the zero 
energy level from the levels just below and above it. 

We next investigate the second equivalence 2 given in 
section III via similar observations. Let us start with 
the honeycomb lattice. The increment in the topological 
number is AI = 1 in the band edge regions, the conven- 
tional Fermi-Landau quantization. This feature indeed 
survives in the adiabatic process t'/t : — * 1 in this re- 
gion. This is the topological equivalence 2. On the other 
hand, for t'/t : — > — 1 these bands merge in pairs to 
produce doubled Chern numbers, which corresponds to 
the fermion doubling due to the multiple band extrema 
discussed in llll CI 

In passing, we mention that one can confirm the 
particle-hole symmetry i.e, the bipartite symmetry in the 
topological numbers, 

I(E F )=-I(-E F ). (15) 



C. From strong to weak magnetic fields 

Having identified the strong-field behavior, let us next 
show how the spectra changes when the magnetic field 
has smaller, more realistic values. In Fig. ^3 we show the 
spectra for the honeycomb lattice. The weaker the exter- 
nal magnetic field becomes, the narrower the gaps as well 
as the band widths become: For <fi — 1 /31, the bands look 
like Landau levels, but the edge states continue to exhibit 
characteristic behavior around the zero energy region (1, 
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FIG. 10: The energy spectra of the honeycomb lattice with zigzag and bearded edges for (a) <j> — 1/5 and (b) (f> — 1/21, (c) is a 
blowup of the (b), and (d) is for <j> = 1/51. The topological number I(Ef) when Ef is in a major energy gaps is also indicated 
next to the vertical axis. 



3, 5, ... edge branches for the consecutive gaps), and 
this region is separated by the van-Hove energies ~ ±1, 
outside of which we recover the free fermionic topologi- 
cal numbers with AI = 1. As discussed above, in strong 
magnetic fields such as <f> — 1/5, only two bands across 
the zero-energy are Dirac-Landau levels, while the next 
bands with energy ~ ±1 are van-Hove bands. Decreas- 
ing magnetic field makes the number of Dirac fermion 
bands increases: The number of edge states in Fig. 1101 is 
precisely in accord with this. 

VI. SUMMARY AND DISCUSSIONS 

In summary, we have shown that the Landau levels are 
divided by the van Hove singularities into two regimes: 
One is effectively described by Dirac particles, and the 
other by ordinary finite-mass fermions. Remarkably, this 
persists as we convert the lattice into square or 7r-flux 
lattices. We have confirmed this both from the Chern 
number argument on Hofstadter's diagram and by the 
bulk-edge correspondence. 

The van Hove singularities have turn out to play a key 
role in separating the Dirac-Landau and Fermi-Landau 
regions. Thus the anomalous QHE is not unique to the 
honeycomb lattice, but continuously shared by a class of 
lattices (except by the square lattice) . The Dirac-Landau 
behavior can also occur when there are multiple extrema 
in the dispersion, as in the region between honeycomb 
and 7r-fiux lattices. 

The properties revealed here for the whole spectrum 
of the honeycomb lattice may be experimentally observ- 
able if the chemical potential can be varied over a wide 
range. In real graphene samples there may be disorder, 
in which case we are talking about a dirty Hofstadter sys- 
tem. However, non-monotonic behaviors should survive 



the disorder as far as the degree of disorder is not too 
large, as has been indicated by a numerical calculation 
for a dirty Hofstadter system—. As for the 7r-flux lattice, 
the model effectively describes the excitation spectrum of 
d-wave superconductors, so there may exist a relevance. 
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APPENDIX A: MOMENTUM 
REPRESENTATION OF THE HAMILTONIAN 

We derive the Hamiltonian for the honeycomb lattice 
in the momentum representation for the flux per hexagon 
<f> = p/q, where p and q are mutually prime integers. We 
label the unit cell as j = j±ei + ]2&2- Let r M be the 
reciprocal vectors (divided by 27r) satisfying e M • r v = 
<5 M „. In an external magnetic field, it is convenient to 
introduce a larger unit cell along coordinate 1 as J = 
Ji(gei) + J2&2- In this unit cell, the fermi operators are 
denoted as c sm (J) with s = »,o and m — 1, ••• ,q. 
Let K be the corresponding momentum defined by K = 
Ki(ri/q) + K<iXi. Then the Fourier transformation can 
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be defined as 



dK 



e lKJ c sm (K). 



Note that £ [0,2w]. Introduce K\jq = k\ and K 2 = 
k 2 , where fc M denotes the momentum in zero magnetic 
field. It is convenient to introduce operators c^(k) — 

(4(fc), d 2 (k), 4(fc), . . . , c t q (k), ci g (k)) . with 

these we obtain the Hamiltonian in the momentum rep- 
resentation Eq. as 



h(k) = 



( d x fx 

A d 2 f 2 



-iqki ft 



FA 



dq-1 f, 
fn-1 



9-1 



where 



( 1 + e -ik 2 +i2^j 

/j( fc 2) = * ( 1 + (t'/t)e-**»+*!hr*y+i/a) ) 



(Al) 



APPENDIX B: GAUGE TRANSFORMATION 
AND THE TRANSFER MATRIX 



In the transfer matrix formalism, the energy bands of 
the bulk system can be determined by the condition 



p 2 - pTiM + dct M = 



(Bl) 



that satisfies \p\ = 1. 

Since the transfer matrix has complex matrix elements, 
let us separate the hopping parameters into the magni- 
tude and the phase as 



where r s is a real non-negative parameter. Obviously, 
r s{ji + 9, k 2 ) — r s (Ji, k 2 ). In what follows, the k 2 depen- 
dences are suppressed for simplicity. 



Then the matrices in Eq. QlOjl can be expressed as 



e i0o.(ii)+i'0.o(ji-i) 

gi0.oO'l-l) 



e i<Po.(jl)+i<P.otil) 



ro.Ui 
1 



r.oQ'i-l) 
ro.(ji) 





e -j0.oO'i-i) 
1 I ' 



E ro.(ji) 

1 



g-#o.O'l) Q 

1 



By multiplying the above two matrices, we can see that the transfer matrix (|llf> can be written as 



gj0o.O'i)+«0.o(ii) 

Q g#o.(jl) 



U(j 1 + l)M t (j 1 )U- l (j 1 ), 



E 



r .(ji) 

r.o(jl) r.o(jl) 



1 







1 







1 







g^«o(jl-l) 



where M t is a real matrix defined by the product of the second and third matrices in the r.h.s. of the first line, and 
U (ji ) is a diagonal matrix defined by 



/ ji-i 

0.0$) 



U{ji) = expi 



Ji 







Here, the sum over j[ starts from some integer, say, —1. fore, we arrive at an expression for the unit cell transfer 
The matrix U(j\) is nothing but a local gauge transfor- 
mation for the wave functions <p {ji ) and ip, (ji ) . There- 
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matrix H3|) . 

M = U(q + 1)MIJ-' L (1), 

where M is a real matrix defined by M = 11^=1 ^t(ji)- 
Note that 

U- 1 (l)U(q + l)=e ie l, 

where is defined by 8 e X^ 1= i S s 'frsih) = 
Argi.o(ji)t .(ji) (mod 2tt). Finally, we have 

detM = e 2ie , 
TrM = e ie TrM, 



where we have used det Mj= 1 in the first line. It should 
be stressed again that TrA/ is real. Thus we conclude 
that the quadratic equation has solutions that satisfy 
\p\ = 1 if |TrM| = | TrM | < 2. 

As a function of the energy, the eigenvalue of the trans- 
fer matrix is written as 

9 = \ ( Tr M ~ ^ Tr M -* 2 ~ 4 det M ] ' 

where we have to fix the branch of the square root. Then 
the Riemann surface, £2,3-1, of the Bloch function is 
given by that of the eigenvalue p{E) and its genus is 
generically 2q — 1 which are the number of energy gaps. 



* The quality of the figures is reduced due to the cond-mat 
file size limit. The original pdf file is obtained from the link, 
http : //pothos . t .u-tokyo . ac . jp/Lab/Research/preprint 
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